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We have used high resolution observations obtained at the AAT with UHRF (R ~ 
"^ I 100,000) and at Geniini-S with b-HROS {R -^ 150,000) to determine magnesium isotope 

V^ I ratios for seven u Cen red giants that cover a range in iron abundance from [Fe/H] = 

^ ; -1.78 to -0.78 dex, and for two red giants in M4 (NGC 6121). The w Cen stars sample 

both the "primordial" (i.e., 0-rich, Na and Al-poor) and the "extreme" (0-depleted, 



Ph! 



2 ' Na and Al-rich) populations in the cluster. The primordial population stars in both 

'^ . oj Cen and M4 show (^^Mg, ^^Mg)/^^Mg isotopic ratios that are consistent with those 

^ found for the primordial population in other globular clusters with similar [Fe/H] values. 

The isotopic ratios for the uj Cen extreme stars are also consistent with those for extreme 

^ I population stars in other clusters. The results for the extreme population stars studied 

rn ' indicate that the ^^Mg/^^Mg ratio is highest at intermediate metallicities ([Fe/H] < 

(y-v ■ -1.4 dex), and for the highest [Al/Fe] values. Further, the relative abundance of ^^Mg 

^^ ■ in the extreme population stars is notably higher than that of ^^Mg, in contrast to 

■^ . model predictions. The ^^Mg/^^Mg isotopic ratio in fact does not show any obvious 

^.f-^ . dependence on either [Fe/H] or [Al/Fe] nor, intriguingly, any obvious difference between 

the primordial and extreme population stars. 

K^ I Subject headings: globular clusters: general; globular clusters: individual (w Centauri, 



S 



NGC 5139; M4, NGC 6121); stars: abundances; stars: Population H 



1. INTRODUCTION 

The stellar system uj Centauri has been known to have an unusual stellar population for almost 
four decades. For example, there is now an extensive body of work which shows that, unlike the sit- 
uation for most globular clusters, the member stars of uj Cen possess a large range in heavy element 



abundance together with distinctive element-to-iron abundance ratios (e.g.. |Johnson &: Pilachowski 



20ld : lMarino et al.ll2012l . and the references therein). Based on photometric and spectroscopic stud- 



ies of lower main s equence stars there is also evide nce for a substantial spread in helium abundance 



in the cluster (e.g.. lNorrisll2004l : lPiotto et al.ll2005l ). Further, spectroscopic studies of stars near the 



main sequence turnoff suggest that uj Cen has an age spread of order 2 Gyr (e.g., iPancino et al 
20 111 , and the references therein). These properties, which are not shared, in toto, by any other 
globular cluster, have led to the suggestion that co Cen is the nucl ea r remnant of a disrupted dwarf 



galaxy that was accreted by the Milky Way (e.g., lFreemanlll993l ). iBekki &: FreemanI (|2003l ) have 
shown that this is dynamically plausible despite the tightly bound and retrograde current orbit of 
the cluster (JDinescu et al.lll999l ). 



In one respect, however, cj Cen is very similar to the other globular clusters of the Milky 
Way and that is the pre sence of the so-called 0-Na anti -correlation arnong the cluster stars 
(jNorris fc Da Costal Il995al ). It is now well-established (e.g.. ICarretta et al.ll2010l ) that essentially 
all Milky Way globular clusters show this phenomenon. It is characterized by the presence of two 
or three populations within a given g lobular cluster: prim ordial, intermediate and, in at least some 



clusters, an extreme population (e.g.. lCarretta et al.ll2009l ). The primordial population shows abun- 
dances of the elements C, N, O, F, Na, Al and Mg relative to Fe that are essentially equivalent to 
those for field halo stars with the same [Fe/H] value. The intermediate population shows depletions 
in O of order 0.3 dex, and depletions in C and F, which are coupled with enhancements in N, Na 
and Al. The extreme population shows even larger depletions in O and also depletions in the Mg 
abundance in addition to the N, Na and Al enhancements. 

The usual interpretation of these abundance correlations is that they are the consequence of 
proton-capture reactions occurrir ig near the H-burning shell i n sorn e ast rophysical source (a s origi- 
nally suggested in this context by lDenisenkov h Denisenkoval (Il990l) andlLaneer et al.l (jl993l )). The 



reactions involved are those of the Ne-Na and Mg-Al cycles (e.g., ICavallo et al.lll996l ). Magnesium 
abundances, and in particular Mg isotope ratios, have played a crucial role in the development of 
theories for the o r igin o f the abundance anomalies. In particular, the outcome of the pioneering 
work of IShetrond (jl996l ) helped rule out individual star "deep r nixing" scenario s as a possible ex- 
planation for the observed abundance differences. Specifically, IShetroneJ (j 19961 ) analysed spectra 
of two 0-rich and four extremely 0-poor red giants in M13. The 0-poor stars were also strongly 
Al-enhanced and notably Mg-depleted. The key result was the recognition that the Mg-depletions 
in the 0-poor stars required depletions in the ^^Mg abundance, with the ^^Mg and ^^Mg isotopes 
having above solar abundance ratios. Since the ^capture reaction on ^^Mg in the Mg-Al cycle 
is strongly temperature dependent, significant ^^Mg depletions occur only for temperatures well 
above those found near the H-burning shell in canonical low mass stellar models. Consequently, 
while evolutionary mixing does occur in red giant branch stars, the primary origin of the abun- 
dance anomalies is now accepted as the result of events occurring during the formation of globular 
clust ers that involve the i nclusion of material pr ocessed in either rap idly rotating massive stars 
(e.g., iDecressin et al.ll2007ll. massive b inary stars (jde Mink et al.ll2009l ) or relatively massive AGB 



st ars fe.g.. IVentura k. p'Antonall2005l ) into a second generation of cluster stars. Subsequent work 
by lYong et all (l2003al. I2OO6II on t he Mg isotope ratios for red giants in the clusters NGC 6752 and 
M13, and by lMelendez &: CohenI (|2009|) for M71, have verified and extended the results of lShetrone 
(|l996l ). 



While there are extensive studies of the 0-Na anti-correlation in a number of clusters with dif- 



ferent [Fe/H] values (e.g.. lCarretta et al.ll2009l ). the presence of a [Fe/H] abundance range in uj Cen 
can be turned to advantage as it allows investigation of the 0-Na anti-correlation phenomenon 
as a function of [Fe/H] in a single (admittedly complex) syste m. The richness of to Cen als o al- 
lows large samples o f stars to be studied. The investigations of I Johnson &: Pilachowskil ([20101) and 
Marino et al.l (|201ll ) show that the 0-Na correlation is present at all [Fe/H] values except perhaps 
for the highest [Fe/H] values ([Fe/H] > -1.0 dex). In particular, the maximum value of [Na/Fe] 
seen in 0-depleted stars increases monotonically with [Fe/H] reaching [Na/Fe] ~ -1-1.0 for the most 
metal-rich objects. In this latter population Na-rich stars are dominant despite the apparent pres- 
ence of a significant range in [0/Fe] values. The situation for [Al/Fe] is somewhat different - the 
maximum value of [Al/Fe] in 0-poor, Na-rich stars is relatively constant at [Al/Fe] ~ -1-1.2 for stars 
with [Fe/H] < -1.3, but for higher metallicities t he maximum value of [Al/Fe] decreases significantly 
to [Al/Fe] ^ -1-0.6 at the highe st metallicities (JNorris &: Da Costal Il995bl : iJohnson &: Pilachowski 



2010 



Da Costa fc Marino 



2011 



Possib l e way s to g enerate these abundanc e patterns are discussed 



in, for example, iJohnson &: Pilachowskil (|2010l ) and lD'Antona et al.l (J201ll ). 

The only extensive study of Mg in to Cen stars, however, is that of lNorris fc Da Costal ( 1995bl 1H 
The results for the stars analy sed by these authors (whose d ata f or Na, Al and O exh ibit the same 
trends with [Fe/H] as those of lJohnson fc Pilachowskil (J2010l ) and lMarino et al.l (J201ll )) reveal only 
a minor population of stars with significant Mg depletions. The stars with Mg-depletions are, as 
expected, strongly en hanced in [Al/Fe] and generally also show significant depletions in [0/Fe]. 
In the terminology of ICarretta et al.l (l2009l^ most belong to the extrem e p opulation. W hile the 



sample is not large compared to those of lJohnson &: Pilachowskil (|2010l ) and lMarino et al.l (|2011 



it is notable that these Mg-depleted stars are all relatively metal-poor: the highest [Fe/H] value is 
-1.4 dex. This is no doub t related to the smaller Al enhancements seen at highe r [Fe/H] values 
(JNorris fc Da Costal [l995bl : I Johnson fc Pilachowskil boiol : ba Costa fc Marinoll201ll ). 



It is clear that further information on the process or processes underlying the 0-Na anti- 
correlation could be obtained by the measurement of Mg isotope ratios in a suitably selected se t 



of u} Cen red giants. For example, in the 6 Mq metal-poor AGB model of I Ventura et al.l (|2011 



Mg-Al cycle burning at the bottom of the convective envelope produces surface Mg abundances 
that are dominated by ^^Mg with low ^^Mg abundances. This is in conflict with the e xisting data 
for Al-rich stars in NGC 6752 and M13 where ^^Mg has a higher abundance than ^^Mg ( Yong et al 



2003al . I2OO6I ) . Confirmation of this discrepancy would likely drive a reassessment of the theoretical 
treatments, such as mixing processes and nuclear reaction rates, employed in the stellar models. 

Against this background, the purpose of the present paper is to determine the Mg isotope ratios 
for a representative set of red giants in w Cen and for two red giants in the "normal" globular cluster 
M4 (NGC 6121). The selection of the sample is described in Section 2, together with a presentation 



^We note for completen ess that ISmith et al.l (2000|) derived Mg abundances for 10 ui Cen r ed gi ants while 
Brown fc WallersteinI (199^) provided Mg abundances for 6 stars, one in common with lSmith et al.l (2000|). 



of the observations and data reduction. The analysis to derive the Mg isotope ratios is discussed 
in Section 3 and the results for the w Cen and M4 stars are presented in Section 4. In Section 5 we 
combine our results with previous studies and discuss the implications in the context of existing 
models for the synthesis of the ^^Mg and ^^Mg isotopes. 



OBSERVATIONS AND REDUCTIONS 



2.1. Sample Selection 

The uj Cen stars to be observed in this progra m were selected from the set of red giant 
members observed at high resolution and analysed bv lNorris fc Da Costal (jl995bl ). These authors 
give atmospheric parameters, iron abundances [Fe/H], and element-to-iron abundance ratios for a 
biased sample of 40 stars that covers nearly the full range of [Fe/H] values seen in the cluster. In 
particular, we first selected ROA 94, which with [Fe/H] = ~l-78 is o r ie of the most meta, l-poor stars 
in w Cen (jNorris fc Da Costal Il995bl : Ijohnson k Pilachowskil boiol : iMarino et ahlboid). The star 
has an enhanced [a/Fe] ratio and low values for [Na/Fe], [Al/Fe] and [Ba/Fe] (|Norris &: Da Costa 
1995br). It is therefor e a member of the primordial population of the cluster in the terminology of 
Carretta et al.l (|2009l ). and with the low [Fe/H] value, a good candidate with which to characterise 
the "initial" value of the Mg isotope ratios in the cluster. The remaining stars in the observing 
list were pairs of stars selected to have similar [Fe/H] values but with significantly different [0/Fe], 
[Na/Fe], etc, abundances. In other words, we selected pairs of primordial and extreme population 
stars with similar [Fe/H] values to investigate the difference in the Mg isotopic ratios between 
the different populations. The set of pairs also encompassed a range of [Fe/H] values so that any 
difference between the populations as a function of [Fe/H] could be explored. Three such pairs of 
stars were subsequently observed: ROA 43 and 100, ROA 132 and 150, and ROA 201 and 248. 
Their abundance characteristics and at mospheric parameters along with those for ROA 94 are given 
in Table [H The values are taken from iNorris &: Da Costal (Il995bl ) . 



We note that the primordial stars ROA 43, ROA 132 and ROA 201 belong to the "CO-strong" 
population in a; Cen, while the extreme stars ROA 150 and ROA 248 belong to the "CO-weak" 
population. No CQ -strei igth classification is availab le for ROA 94 or ROA 100. As discussed in 
Persson et al.l (jl980l ) and lNorris &: Da Costal (jl995bl . Section 2.1), the w Cen CO-strong stars have 
stronger near- I R CO indices than those of all of the red giants in other globular clusters for which 
Persson et al.l (Il980l ) present data, while the CO-weak stars have CO indices comparable to those 
of other globular cluster red giants at similar [Fe/H] values. For com pleteness we note a l so that 
the extreme stars ROA 100 and ROA 150 are classified "CN-strong" in lNorris k. Da Costal (Il995bl ) 
while the primordial star ROA 94 is classified "CN-weak". No CN-band strength classification 
information is available for the primordial stars ROA 43, ROA 132 and ROA 201 and for the 
extreme star ROA 248. 



Six of our seven stars have also been observed in I Johnson k Pilachowskil (J2010l ): stars ROA 94 



43, 100, 132, 150 and 248 correspond to Ijohnson fc Pilachowskil (|20ld ) stars 39037, 10012, 44449, 
55114, 73025 and 55071 , respectively. For these six st ars we find mean [Fe/H] , [Q/F e] and [Na/Fe] 



differences, in the sense iNorris &: Da Costal (|l995bl ) - I Johnson &: Pilachowskil (J2010l ) of 0.00, -0.01 



and +0.03 dex, respectively. The standard deviations are 0.10, 0.13 and 0.09 de x, which are consis- 



tent w ith the combined uncertainties in the abundance determinations. Although lJohnson &: Pilachowski 



201C ) do not list [Al/Fe] ratios for the six common stars, [Al/Fe] values are giv en in 



20081). Using those data we find a mean difference in [Al/Fe] in the sense iNorris &: Da Costa 



Johnson et al. 



1995bl ) - iJohnson et al.l (J2008l ) of -0.08 dex thou gh with a disturbin g ly large standard deviation 



of 0.47 dex. Since the [Al/Fe] determ inations of INorris &: Da Costal (jl995bl ) are in good agree- 
ment with those in ISmith et al.l (J200d ) (mean difference of -0.07 dex with a = 0.16 dex for 3 
stars which co ver an [Al/Fe] r ange of ~1 dex), we ascribe the large a primarily to large uncer- 



tainties in the 



Johnson et al. 



that the Johnson et al. 



tinct from R ~ 18,000 for 



(|2008l ) [Al/Fe] determinations. In that context it is worth noting 
(l2008l) spectra have only modera te spectral resolution R ^ 13,000 as dis- 



Johnson fc Pilachowskil (|2010l ) and R ~ 38,000 for INorris fc Da Costa 



(Il995bl). Sirn i larly, star ROA 132 is the only one from our sample in common with the work of 



Marino et al 



(I2OIII I ; it cor re sponds to the i r star 206284. The abundance differences in the sense 



Norris fc Da Costal (Il995bl l - iMarino et aD (|201lh are 0.19, -0.18 and +0.07 dex for [Fe/H], [0/Fe] 
and [Na/Fe] , respectively. Unfor tunately, this star is not among those with C, N and O abundance 



measures m 



Marino et all (|2012l ). 



The two M4 stars were selected as the brightest pair of "CN-weak" (L151 4) and "CN- s trong " 
(L2406) stars in the detailed abundance study of M4 red giant s carri ed out by llvans et al.l (|l999l ). 
Based on the [0/Fe] and [Na/Fe] values given in llvans et al.l (119991 ) . star L1514 is a member of 
the primordial population of M4 while L2406, which has an oxygen abundance 0.22 dex lower 
than L1514 and a sodium abundance 0.30 dex higher, is a member of the cluster's intermediate 
population. The two stars differ marginally in [Al/Fe] with L2614 0.11 dex highe r, but have 



essent ially identical Mg abundances. Indeed in the sample of M4 red giants observed by llvans et al 



(j 19991 ). there is no evidence for any intrinsic range in [Mg/Fe] or any obvious correlation between 



for the two stars from 



Ivans et al 



[Mg/Fe] and [Na/Fe] or [Al/Fe] (Ivans et al.lll999l ). The atmospheric parameters and abundances 



(|l999l ) are given in Table [TJ T he two M4 stars we re also included 



in the stud y of lYong et al.l (|2008al ) who found similar results to llvans et al.l (jl999l ). In particular, 
Yong et al.l ()2008al ) find an oxygen abundance difference of 0.18 dex (in the sense L1514 - L2406) 
and [Na/Fe], [Mg/Fe] and [Al/Fe] differences of -0.36, -0.08 and -0.03 dex, respectively. The 
oxygen and sodium abundances confirm the population classification for these stars. 



Star L1514 is also included in th e M4 study of 



the oxygen a nd sodium abundan ces of 



Marino et al 



Marino et al.l ((20081), but L2406 is not. Using 



(I2OO8I ) , L1514 is again class ified as a primordial 



star. As did llvans et al.l (|l999l ) and I Yong et all (|2008al ). iMarino et al.l (|2008l ) find no significant 
magnesium - aluminium anti-correlation, and that [Mg/Fe] does not correlate with either [0/Fe] 
or (anti) correlate with [Na/Fe] within t heir extensive samp le of M4 red giant members. The total 
[Al/Fe] range is also small at <0.4 dex (IMarino et al.ll2008l ). 



Table 1. Globular Cluster Program Stars 



Object Teff log g [Fe/H] [C/Fe] [0/Fe] [Na/Fe] [Mg/Fe] [Al/Fe] Pop Type 



w Cen ROA 94 4200 0.7 -1.78 



-0.67 0.25 



-0.26 



0.51 <0.17 



u Cen ROA 43 3950 0.4 
CO Cen ROA 100 4150 0.7 



-1.47 -0.21 0.39 0.21 0.52 0.06 P 

-1.49 <-0.54 <-0.43 0.58 0.05 1.15 E 



u Cen ROA 132 3900 0.3 -1.37 -0.22 0.35 0.17 0.56 -0.32 P 

u Cen ROA 150 3950 0.6 -1.25 -0.83 -0.46 0.59 0.26 1.06 E 



u Cen ROA 201 3750 0.5 -0.85 
u Cen ROA 248 3850 0.6 -0.78 



-0.30 



0.37 
-0.33 



-0.12 
0.75 



0.45 
0.24 



-0.07 
0.68 



P 

E 



NGC 6121 L1514 3875 0.3 -1.18 -1.00 0.41 0.01 0.40 0.44 P 

NGC 6121 L2406 4100 0.5 -1.18 -0.84 0.19 0.31 0.37 0.55 I 



NGC 6752 mgO 3928 0.3 



-1.62 



-0.08 



0.65 



0.46 



1.06 



E 



Note . — Except for the last column, the data for the to Cen entr ies are taken from lNorris fc Da Costa 



(|l995bl ). while those f or the M4 stars come from llvans et al.l ()1999l ). The NGC 6752 mgO entries use the 
results of lYong et al.l (J2003al ). including the solar abundances given in that paper. In the last column "P" 
indicates the star is classed as belonging to the primordial population, "I" for the inte r media te population 
and "E" for stars in the extreme population, using the terminology of ICarretta et al.l (|2009l ). uj Cen stars 
ROA 132 and ROA 150 were observed with UHRF at the AAT as were the two M4 (NGC 6121) stars. 
UJ Cen stars ROA 94, 100, 201 and 248, together with NGC 6752 mgO, were observed with b-HROS at 
Gemini-S. Star ROA 43 was observed at both facilities. 



2.2. Observations 

The observations were obtained with two different telescope/spectrograph combinations. The 
first employed the "Ultra-High Resolution Facility" (UHRF jj on the Anglo- Australian Telescope 
(AAT). The instrument was used in the resolution 300,000 mode but with a 0.6" slit in place of 
the normal imager sheer. This setup yields spectra with an effective resolution R ~ 100,000 and 
a scale of 0.018A per binned pixel with the Tek Ikxlk detector in use at the time. Due to the 
limited wavelength coverage available with this detector, the central wavelength of the spectrograph 
setup had to be adjusted for the geocentric velocity of each target to ensure the isotopic abundance 
sensitive MgH features (A ~ 5134.55 - 5140. 2A) were centred. Data were obtained in two runs in 
2000 April and 2001 March under seeing conditions that varied from 1.2" to 2.4", with a typical 
value of 1.8". As the small slit size compared to the seeing limits the instrument efficiency, total 
exposure times of many hours (broken into individual 1800 sec integrations) were needed to achieve 
signal-to-noise (S/N) values that would be useful for this project. In the event, as we shall see 
below, we were able to obtain S/N values of ~50-60 per pixel for the a; Cen stars ROA 43, 132 
and 150 with this instrument setup for total exposure time of 10, 10 and 17.5 hrs, respectively. For 
the M4 stars L1514 and L2406, the S/N values were ~70 and ~45 for total exposure times of 2.0 
and 3.5 hours, respectively. We also obtained s pectr a of bright field stars for which Mg isotope 



ratios have been determined by lGay fc LambertI (J2000l ). Flat field, Th-Ar arc lamp and dark frame 



exposures were obtained each night during both observing runs. 

The second telescope/instrument combination used was the (now decommissioned) b-HR0q3, 
bench-mounted fiber-fed high dispersion spectrograph on the Gemini-S 8m telescope. The spec- 
trograph was located in the pier of the telescope and fed via fibres from a cassette mounted in 
GMOS-S that locates them in the focal plane of the telescope. Guiding and tip/tilt image motion 
compensation is achieved through use of the GMOS-S On-Instrument Wavefront Sensor. We used 
the 'object only' mode in which a single fibre with a projected diameter of 0.9" on the sky transfers 
the star light to the spectrograph. An image-slicer at the spectrograph entrance slit then reformats 
the fibre output into seven 0.14" wide slices with a nominal resolution of i? ~ 150,000. The central 
wavelength of the echelle was set to A5135A in order 40. With the single 2048 x 4603 pixel (spatial 
X wavelength) CCD detector this provided sufficient wavelength coverage that no adjustment in 
the central wavelength was required to compensate for the range of radial velocities exhibited by 
the various targets. On-chip binning of x4 in the spatial direction was used and the scale in the 
spectral direction was O.OIOA per pixel. 

The observations were carried out during a two night (2006 UT March 16 and 18) classical 
run. Both nights were clear with excellent seeing (~0.6"). Typical total exposures for the program 
targets were of order 2 hours per star, broken into individual 1800 sec integrations. Observations 

^ www. aao.gov.au/astro/uhrf. html 

■^ www. gemini.edu/sciops/instruments/hros/hrosindex. html, accessible via www.gemini.edu/node/10002 



were obtained of lo Cen red giants ROA 43, 94, 100, 201 and 248, where we note that ROA 43 
was also observed with the AAT/UHRF set up. In ad dition, we also obse rved the globular cluster 



red giant NGC 6752 mgO which has been analysed bv lYong et alj (|2003al ). whose data we present 



for comparison purposes in the final row of Table [TJ Spectra wer e also obtained o f three bright 



field stars f or whi ch Mg isotope ratios have been determined by iGav &: LambertI (|2000l ) or by 



Yong et alJ (l2003al ). Flat field and arc lamp exposures were also obtained each night using the G- 
CAL facility calibration unit. All frames were bias-subtracted and trimmed of the overscan regions 
at the telescope. 



2.3. Data Reduction 

The data for the AAT/UHRF and Gemini-S/b-HROS configurations were reduced by using 
the FIGARO packagqj- Individual spectra were flatfielded and wavelength calibrated with spectra 
of the quartz and arc lamps, respectively; they were cosmic-ray cleaned using "bclean"; curvature 
in the spectra was removed with "sdist" and "cdist" ; individual observations of each object were 
wavelength shifted to bring them into the restframe; and for each object the individual spectra were 
co-added, sky-subtracted, and finally continuum normalized by using spectra of a rapidly rotating 
B star (/3 Centauri) to constrain the shape of the instrumental response function. An essential 
difference between the two instrumental configurations is that for UHRF the spectral lines were 
set perpendicular to the dispersion direction on the detector, while for b-HROS this was not the 
case. For the latter, we used curvature-corrected arc spectra to determine and remove the slope 
of the spectral lines in the spatial direction by estimating the differential wavelength shifts in that 
direction, and shifting the stellar spectra by the appropriate amounts before summing the counts. 

For the three oj Cen stars observed at the AAT with UHRF the net counts in the continuum 
near 5140A he in the range 2700-3100 electrons per O.OISA pixel. For M4 L1514 and L2406, the 
net counts were 4800 and 2200 electrons per pixel, respectively. This suggests the S/N values for 
these AAT globular cluster red giant spectra lie in the range ~50-70 per pixel. For the Gemini-S 
data the net counts range from 6300 (ROA 43) to 15900 (NGC 6752-mgO) electrons per O.OIOA 
pixel, suggesting S/N in the range ~80-130 per pixel. These values are confirmed by a direct 
determination of the RMS dispersion of the counts in the stellar continuum in small bands of width 
'^0.3A at 5135. 5A and 5137. 9A in the spectra of the hotter and higher S/N uj Cen stars ROA 94 
and 100 and in NGC 6752 mgO. The RMS dispersions lead to S/N values in the range 90-130, 
consistent with those calculated from the net counts. 

Figures [J and [2] show the final reduced spectra for the uj Cen stars over wavelength regions that 
include electronic transitions of the MgH molecule, and the features specific to each Mg isotope are 
identified. Similarly, Fig. [3]shows the isotopic ratio sensitive regions for the two M4 stars. Apparent 



'www.aao.gov.au/figaro 



differences in the isotopic compositions are immediately evident from these figures, especially for 
the uj Cen stars. We now seek to quantify the differences. 



3. ABUNDANCE ANALYSIS 



The isotopic abundances o f magnesium were measured using the identical approach to that in 
Yong et al.l (|2003al . |2004| . |2006| ) and thus can be briefly describ ed here. We u sed one dimensional 
local thermodynamic equilibrium (L TE) model atmospheres frq m lKuruca (119961) and ado pted stellar 
parameters for the u: Cen s t ars fro m lNorris &: Da Costal (jl995bl ). from llvans et al.l (j 19991 ) for the M4 
stars and from lYqng et al.l (l2003al) for NGC 6752 mgO. The particular stellar parameters employed 



are not crucial as I Yong et al.l (J2004l ) has shown that the derived isotope ratios appear insensitive to 
the adopted stellar parameters, non-LTE effects and inadequacies in the model atmospheres. We 
an alysed the three lines from th e A — X electronic transition of the MgH molecule recommended 
bv JMc William &: Lambert! (|l988l ) as providing the most reliable isotope ratios (the numerous other 
MgH lines present in the spectra are more strongly affected by known, and unknown, blends). The 
three lines are (i) 5134. 6A from the (5i(23) and i?2(ll) lines of the (0,0) band, {ii) 5138. 7A from 
the (0,0) Qi(22) and (1,1) Q2(14) lines, and (m) 5140.2A from the (0,0) -Ri(lO) and (1,1) i?2(4) 
lines. 



We note that other authors choose to analyse different MgH features - for example. JMelendez &: Cohen 



( 20091 ) prefer the 5134. 3A feature to that at 5134. 6A - potentially resulting in systematically differ- 
ent isotopic ratios. We therefore make no strong claim for the validity of the absolute values of our 
isotopic ratios but argue that consistent analysis allows meaningful comparisons in a relative sense. 
It is also likely that known and unknown blends affecting the MgH analysis are more problematic 
for cooler and more metal-rich stars, potentially leading to feature-to-feature differences in the 
derived isotopic ratios. Non-LTE and 3D effects may also become more important for cooler tem- 
peratures. Some care must then be taken in comparing results as a function of metallicity although 
as shown below, our results are consistent with published values. For example, our Mg isotopic 
abundance determination for the cool red giant HR 2140 (Tpfj = 4000 K, log g = 0.9, [ Fe/H] = -0.9 



(jMcWilliam k, Lambertlll988l ) ) is consistent with that of JMc William fc LambertI (jl988l ) . Indeed the 
comparison of our results with published values that use the same MgH wavelength regions argues 
that our results are internally consistent, although a relative offset cannot be ruled out. Neverthe- 
less, additional studies seeking to identify blends affecting these and other MgH features would be 
welcome. 



Synthetic spectra were generated using the LTE stellar line analysis program MOOG (|Sneden 



19731 ). with the line broadening, assumed to have a Gaussian form, estimated for each star by fitting 
the profile of the 5115. 4A Nil line. The values of the line broadening employed range from 6 to 
10 km s~^, with an uncertainty of order 1-2 km s~^. Each of the three MgH lines was analysed 
independently with three free parameters: the ^^Mg/^^Mg and ^^Mg/^^Mg isotopic ratios and the 
total Mg abundance, loge(Mg). A large range of parameter space was explored and the optimum 
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Fig. 1. — Spectra for the uj Cen program stars for wavelengths between 5133A and 5136A. The 
spectra are ordered approximately by metallicity with ROA 94 being the most metal-poor and 
ROA 248 the most metal-rich, respectively. The primordial population star spectra are plotted as 
red dots, while those for the extreme population stars are plotted as black lines. The positions of 
the ^^MgH, ^^MgH, and ^^MgH features are marked by vertical lines. 
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Fig. 2. — As for Fig. [T] except that the wavelength region shown is 5138A to 5141A. 
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Fig. 3. — The upper panel shows the wavelength region 5133A to 5136A and the lower panel shows 
the wavelength region 5138A to 5141A for the M4 primordial population star L1514 (red dots) 
and for the M4 intermediate population star L2604 (black dashed line). As for Figs. [1] and [2] the 
positions of the ^^MgH, ^^MgH, and ^^MgH lines are marked by vertical lines. 



13 



value for each parameter identified by locating the x^ minima. Following iBevington &: Robinson 
(jl992l ). the 1-cr confidence limit for a given parameter was taken as the value at which A^^ = x^ 
— Xmin ~ -*-• Thus, for each line in each star, we obtain an optimum value and an uncertainty for 
25Mg/24Mg and ^SMg/^^Mg. The final ratio ^^Mgr^SMgi^^Mg was then calculated as the weighted 
mean of the three determinations. The individual determinations and the final adopted values for 
the uj Cen and M4 stars, and NGC 6752 mgO, are given in Table O 

When expressing the isotopic ratios as ^^Mg:^^Mg:^^Mg = (100 — b — c):b:c, the formal errors 
in the determinations from each region are small; 6 it 1 and c it 1. However, inspection of the values 
in Table [21 for example, suggests that the real uncertainties are larger, being driven by systematic 
uncertainties in the continuum fitting, the adopted macroturbulent velocities, and (known and 
unknown) blends. As discussed in previous studies of stellar Mg isotope ratios, the ^^Mg/^ Mg 
ratio is often more uncertain than the ^^Mg/^^Mg ratio since the isotopic shift is smaller, and as 
a result, the ^^MgH feature is less well separated from the (generally) dominant ^^MgH line. This 
effect can be seen qualitatively in what follows in Figures El El and [7| and quantitatively from the 
X^ analysis. We also note that the isotope ratios derived from each of the three regions are usually, 
but not always, in good agreement. In some rare cases (e.g., u Cen ROA 201), the ratios from the 
three regions show poor agreement for which we have no obvious explanation. 

An additional assessment of the uncertainties can be obtained from independent analysis of 
repeat observations, and from a comparison of isotopic ratios found here with the results of oth- 
ers. For example, three separate spectra of Arcturus, obtained during the AAT/UHRF runs, give 
24]y[g.25]y[g.26]Yjg ratio s of 81:9:10, 84:6:10 and 80:10 :10 which are ah in good agreement with the 
value 82:9:9 found bv iMcWilliam k Lambert! (|l988l ). There are also three separate AAT/UHRF 
spectra of the star HR 2140. As for Arcturus, our determinations for this star of 80:10:10, 
82:8:10, and 83:7:10 are in excellent accord with each other and with the value 88:6:6 given in 



McWilliam &: LambertI (119881 ). Similarly, during the Gemini-S/b-HROS run three field stars with 
previous Mg isotope ratio determinations were observed. The stars are HD 26575 (HR 1299), 
HD 64606 and HD 103036 and we find isotopic ratios of 81:5:14, 89:2:9 and 87:2:11, respectively. 
The published isotopi c ratio determ inations for the firs t two are 72:18:10 and 82:10:8, respectively, 
from lOav fc Lambert! (|200nl ). while lYong et al.l (|2003al ) hst 94:0:6 for HD 103036. Our values are 
consistent with the previous determinations wit h the possible ex c eption of the ^^Mg/^'^Mg ratio 
for ~solar abundance star HD 26575 where the iGay Sz LambertI (l200d ) valuqj is notably higher 
than the value from our analysis. We note, howe ver, that our isotopic ratios f or this star are 
consistent with the 80:10:10 determination given bv lMc William &: LambertI (|l988l ). Perhaps more 
significantly, we find from our Gemini-S/b-HROS observations an isotope ratio of 62: 9:29 for the 



star NGG 6752 mgO, which is in excellent agreement with the value of 63:8:30 given in lYong et al 



^The value o f ^'^Ms; /^'*MK for HD 26575 (0.25) is not shown in the the ^^Mg/^^Mg versus [Fe/H] diagram of 
Gay fc LambertI (2000|) (their Figure 8) as it hes outside the range depicted, notably higher than the ^^Mg/'^^Mg 



values for stars of comparable metallicity. 
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(J2003al ) from VLT/UVES spectra. The isotopic ratios derived from the three wavelength regions 
and the mean values for this star are also given in Table [2j Overall, based on these comparisons, 
we estimate that the uncertainties in our determinations are at the level 6 it 4 and c =b 4 when 
writing the ratio as ^^Mgi^SMgi^^Mg = (100 - b - c):b:c. 



RESULTS 



In the followi ng we will r nake use of the results of lYong et alj (J2003al ) for the red giants in 
NGC 6752 and of lYong et all (|2006l ) for M13 and M71, which have been derived u sing essentially 



the sa me procedure as that adopted here. For NGC 6752, using the terminolog y of ICarretta et al. 



(J2009l ) and the sodium and oxygen abundances given in lYong et al.l (J2003al ). stars NGC 6752 
mgl2 and NGC 6752 mg24 are classified as members of the primordial population, while stars 
NGC 6752 mgO, mg21, mg22 and 702 are members of the extreme population. The other 14 
red giants are members of the intermediate population. Similarly, for the four M13 red giants 
studied in lYong et al.l (J2006l ) , and again using the sodium and oxygen abundances from that work, 
star L598 is a member of the primordial population, star L629 is a member of the intermediate 
population and stars L70 and L973 both b elong to the extre me population. These 4 stars are a 
subset of those studied by lShetrond ( 19961 ). lYong et al.l ( 20061 ) separate ^^Mg from ^^Mg but both 
studies find si milar (^^Mg+^^M g)/^'^Mg ratios for a given star. Further, the si ngle M71 star also 



analysed in lYong et al.l (1200611 is a CN-weak object (IMelendez &: CohenI |2009| ) with sodium and 
oxygen abundances (JYong et al.ll2006l : IMelendez &: Cohenll2009l ) that classify it as a member of the 
primordial population of the cluster. The sodium-to-iron, magnesium-to-iron and aluminium-to- 
iron abundance ratios as a function of [0/Fe] for all these stars are shown in Fig. HI along with the 
equivalent data for the uj Cen and M4 stars from Table [Ij 



4.1. UJ Cen ROA 94 

For this metal-poor star we find ^^Mg:^^Mg:^^Mg values of 88:2:9 with the standard uncertainty 
in the ratios of 6±4 and c±4 where ^^Mg:^^Mg:^^Mg = (100 — 6 — c) : 6 : c. The relative abundances 
of ^^Mg and ^^Mg found here for this star are somewhat higher, partic ularly for ^^Mg, than thos e 



for field halo dwarfs at this metallicity, which are typically 97:1:2 (e.g.. IMelendez &: CohenI 120071 ) . 
However, the relative Mg isotope abundances appear to be some what lower than th e ~80:10:10 



values found for the primordial red giants in NGC 6752 and M13 (JYong et al.ll2003al . l2006l ) . The 
[Fe/H] values for these clusters are higher, by ~0.2 dex, than that of ROA 94. 
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Fig. 4. — Sodium-to-iron (upper panel), magnesium-to-iron (middle panel), and aluminium-to- 
iron (lower panel) abundance ratios are plotted against the oxygen-to-iron abundance ratio. The 
three extreme population stars in co Cen are shown as red filled circles while the four primor- 
dial population red giants are plotted as red open circles. The abundance ratios are taken from 
Norris fc Da Costal (jl995bl ). The M4 stars L1514 (primordial population, blue square) and L240 6 
(intermediate populatio n, blue six-pointed star) are plotted using the data of llvans et al.l (|l999l ) . 
The NGC 6752 data of lYong et al.l (J2003al ) are plotted as filled stars (extreme population), open 
5-point stars ( primordial popula tion) and plus-signs (intermediate population). Similarly, using 
the data from lYong et al.l (120061 ) . the two extreme population stars in M13 are plotted as filled 



triangles, the single primordial population star is shown as an open tria ngle and the inter mediate 
population star is plotted as an x-symbol. The sole M71 star analysed in lYong et al.l ( 20061 ). which 
is a primordial population star, is plotted as an open diamond. 
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4.2. uj Cen ROA 43/100 and ROA 132/150 

As the derived Mg-isotope ratios for these two pairs of stars are found to be similar (see Table 
[2]), we discuss them together. Figures [5] and [6] show our fits to the 5134. 6A MgH feature; the 
fits to the 5138. 7A and 5140. 2A MgH features are similar. For the primordial stars ROA 43 and 
ROA 132 we find essentially identical relative abundances of ^^Mgr^^Mgt^^Mg = 76:10:14, with 
again uncertainties of order ±4 in the b and c values. These values are in reasonable agreement 



with the ^80 :10:10 values found for the primordial stars in NGC 6752 and M13 by lYong et al. 



(|2003al . I2OO6I ). This is perhaps not surprising given that the [Fe/H] values for NGC 6752 and M13 
are only ~0.1 to 0.2 dex lower than the [Fe/H] values of the oj Cen stars, which are given in Table 

m 

For the two extreme population stars we find 2'^Mg:^^Mg:2^Mg = 48:13:39 for ROA 100 and 
24]y[g.25]y[g.26]y[g _ 54; 17:30 for ROA 150. Here, given the larger contribution to the features from 
the heavier isotopes and the lower overall Mg abundance (see Tabled]), the uncertainties in the 
b and c values are somewhat larger, of order ±10. Nevertheless, the requirement for significant 
contributions from the heavier isotopes to permit a good fit to the observations is readily apparent 
in the figures. The isotope ratios for these uj Cen extreme population sta rs are also similar to thos e 



found for the extreme populations in the clusters NGC 6752 and M13 (|Yong et al.ll2003al . l2006l ). 
namely 62:7:31 for NGC 6752 (average of 4 stars) and 51:14:35 for M13 (average of 2 stars). Again 
the iron abundances for these two clusters are only ~0. 1-0.3 dex lower than those of ROA 100 and 
ROA 150 (see Table!]). 

In both pairs the primary difference between the primordial and extreme stars is a decrease in 
the ^"^Mg relative abundance and an increase in the ^^Mg relative abundance; the increase in the 
25]\/[g relative abundance appears more modest. We note also that the difference in total magnesium 
abundance between the primordial and extreme stars in each pair, derived as part of the fitting 
pr ocess, is consistent at the 0.1-0.15 dex level with the differences inferred from Tabled) As argued 



by lYong et al.l (|2004l ) the total Mg abundances derived from the fits to the MgH features are not 
as reliable as those determined from individual Mgi lines; nevertheless, the overall agreement is 
reassuring. 



4.3. u Cen ROA 201/248 

For the 0-rich, Na-poor primordial population star uj Cen ROA 201 ([Fe/H] = -0.85) our anal- 
ysis yields ^^Mg:^^Mg:^®Mg = 74:19:7 while for the 0-poor, Na-rich extreme population star uj Cen 
ROA 248 ([Fe/H] = -0.78) the corresponding values are 77:10:13. Given that the uncertainties 
in the derived ratios are somewhat larger for these cool relatively metal-rich stars, it is not clear 
whether we are seeing any significant difference in isotopic ratios. The synthetic spectrum fits to 
the 5134. 6A MgH feature in the observed spectra are shown in Fig. [7| 



17 





_ 




1 ■ : ■ 






1.00 












/ 


/«"«^4**^ 




- •* 


/ 


\ 


/ 


/# ^'^ 


' / 


/jr ^s5»*^* • 


0.95 


r ^ 


/:■ 


'A 


/ 
/ 




;i 


" 


0.90 


^ '\ 




\ ' 


'Itl 


w * \ 

/ * \ 
• \ 


/ 




ROA 43 J 


0.85 


'- \ 






I 


• 


/ 


« 


79-19-1R = 




: A 


' 


1 'i/ 


• 






100:0:0 I 


^ 0.80 


'_ •< 




1 ' '^ 


• 


• 




7R:q-1.3 


1>\ 

-1—' 


- 




1 '/ 




\.- 




Rfi;15;iq 


£ 0.75 


7 
















-t—' 


- 














- 


c 


















- 0.70 


- 




1 ill 




1 1 1 






1 


(U 


- 




1 1 




' ' 1 




! 


• _ ~ 


■^ 1.00 


5rx 




" V. 


^ 


" _y^» ^\ 




•" 


T^^P^fsi^^^v^—'^ 


o 


^ *«\ 


/> ' 


- jiA •. w 


/ 


jr^ •• \ 


/ 




m • * • *• - 


(U 


: *\ 


!»%«! 


jr\y > 






«^ 


^ 


' •*.>. . 


0^ 0.95 




"^^ 


^ \.:i» 


— ^ * 


# • • 




• 


- 




- 




X?^3i 


/ * 




-y 


■ 


0.90 


|- 






'^ 


/ " 


• 
% 


• 


ROA 100 J 


0.85 


- 














49.1fi-49 ■ 
100:0:0 : 


0.80 


— 








26 




1 


4R-i.^-.^q - 




.^fi-iq-45 ■ 


0.75 


- 




125 




' 








' 


124 


1 


- 


0.70 






1 


1 . . 


1 . . . " 



5134.0 



5134.5 



5135.0 
Wavelength(A) 



5135.5 



Fig. 5. — Spectra of w Cen ROA 43 (upper panel) and ROA 100 (lower panel) compared with 
synthetic spectra. The observations are shown as filled circles and synthetic spectra with different 
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^^Mg mixture, and the thin colored lines show unsatisfactory ratios. The locations of the ^^MgH, 
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Fig. 6. — As for Figure [5] except for oj Cen stars ROA 132 (upper panel) and ROA 150 (lower 
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The lack of obvious diffe rence in Mg isotopi c comp osition between these two oj Cen stars is 
consistent with the results of lMelendez &: CohenI (|2009l ) who found only a small difference in Mg 
isotopic ratios between the CN-weak (primordial) and CN-strong (intermediate) populations in the 
relatively metal-rich cluster M71. For example, the ^^Mg/^^ Mg isotopic ratio for the CN-strong 
stars exceeds that for the CN-weak stars by only ~4 percent (JMelendez &: CohenI |2009| ). The two 
M71 populations in their sample also do not show any significant difference as regards total Mg 
abundance, and the ranges in O, Na and Al abunda nces are small although t hey do exhibit the 
expected 0-Na anti-correlation a nd Na-Al correlation ([Melendez &: Cohenll2009l ) . The [Fe/H] value 
for M71 f-0.91.lYong et al.ll2006l ). which we adopt here, is very similar to those for the u Cen stars. 
Melendez k CohenI (|20oi ) give [Fe/H] = -0.80 for M71. 



4.4. M4 L1514/L2406 

Figure [8] shows the spectrum synthesis fits to the 5134. 6A MgH feature for the two M4 stars. 
The fits for the three wavelength regions yield mean isotopic ratios of 82:6:12 for the primordial 
star L1514 and 80:6:15 for the intermediate population star L2406. Given the uncertainties in the 
determinations there is no apparent difference in the derived Mg isotope ratios. lYong et al.l (l2008bl ) 
derived preliminary Mg isotope ratios for a n umber of M4 r ed giants. Because their spectra had 
a resolution of ~60,000, they assumed (as did lShetronelll996l ) that ^^Mg = ^^Mg. For L1514 they 
found an isotope ratio of 80: 10:10, which is cons istent with our result. An isotope ratio was not 
derived for L2406 but overall lYong et al.l (|2008bl ) did not find any significant evidence for a range 
in the Mg isotope ratios of their M4 stars. 

The lack of any significant Mg isotope ratio difference for the two M4 stars observed here 
is perhaps not surprising given the l ack of any differe r ice in [Mg/Fe] and the (at most) small 
difference in [Al/Fe] for these stars (llvans et al.l Il999l : lYong et al.l l2008al ). Again this is con- 
sistent with the emerging general picture that while the 0-Na anti-correlation is present at all 
metallicities, the extent of Al and esp e cially Mg variati o ns are much reduc e d at higher metal 
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DISCUSSION 



As discussed in iKobayashi et al. 



(2011 



see also IWoosley &: Weaveii ()1995l ): iTimmes et al. 
(|l995l )). the heavy isotopes of magnesium have their origin in a number of nucleosynthetic pro- 
cesses. For example, explosive nucleosynthesis in massive star supernovae can produce ^^Mg and 
^^Mg through the He-burning reactions ^^Ne(a, n)^^Mg and ^^Ne(Q;, 7)^^Mg. Being "secondary" 
isotopes, the production of the neutron-rich isotopes increases with the metallicity of the super- 



20 



c 



> 

-I—' 

o 

01 




5134.0 



5134.5 



5135.0 
Wavelength(A) 



5135.5 
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novae, although so me primary production also occurs (JKobavashi et al.ll201ll ). The nucleosynthetic 
model predictions (JKobavashi et al.ll201ll ) are at least qualitatively consistent with the observa- 
tions of the variati on of the Mg isotopic abundanc e ratios with [Fe/H] in the field star samples of 
Yong et all (|2003bh and lMelendez fc Cohen! (|2007l l. 



Relatively massive (M > 4 M.p)) AGB stars can al s o contribute to the production of the h eavy 



magnesium isotopes (e.g., iKarakas &: Lattanzid l2003l : iKarakasI l20ld : iKobavashi et al.l 12011 



In 

such stars there are a number of production sites: (1) the vicinity of the He-burning shell via 
a-capture reactions on ^^Ne; (2) via the Mg-Al cycle in the vicinity of the H-burning shell; and 
(3) at the b ase of the convective envel ope during the hot-bottom burning phase, again via the 
Mg-Al cycle (JKarakas &: Lattanzidl2003l ). In all cases the processed material can be brought to the 
surface layers and thus can contribute to the e nrichment of the surroundin g inte rstellar mediuin 
via mass loss in a stellar wind. As discussed in iKarakas Sz Lattanzid (120031 ) and iKarakasI ( 2010 ) 
the yields of ^^Mg and ^^Mg are complex functions of the initial abundance of the AGB star and 
of the parameters that control the thermal pulse structure, the extent of thermal dredge-up, the 
conditions in the convective envelope and the rate of mass loss. In the context of the current work it 
is worth noting that the production of ^^Mg and ^^Mg via a-capture reactions on ^^Ne is unrelated 
to the abundance of aluminium. Consequently, since we find that the relative ^^Mg abundance is 
tightly coupled to that of Al, at least for stars with [Al/Fe] > 0.5 (see top panel of Fig. [TO|) . we can 
conclude that it is unlikely that the shell He-burning production site is a significant contributor. 
Indeed it is likely that the primary nucleosynthesis site is hot-bottom burning in relatively massive 
AGB stars. 



In the panels of Fig. [9] we plot the individual isotopic fract ions relative to the t otal M g abun- 
dance against [Fe/H] for the stars observed in oj Cen, using the iNorris &: Da Costal (Il995d ) [Fe/H] 
values, and for the M4 stars using the cluster mean [Fe/H] fron illvans et al.l (119991 1. Also shown 



are t he red giants obse rved in the globular cluster s NGC 6752 (lYong et al. 



M71 (JYong et al.l 120061 ). Since the observations in I Yong et al.l (J2003al . l200d ) were analysed using 



2003a|) and M13 and 



very similar techniques to those employed here, no systematic effects are expected and the data can 
be combineqj. In the figure, filled symbols represent extreme population stars while open symbols 
represent primordial populations stars. Fig. [10] shows the same relative abundance fractions but in 
this case plotted against the [Al/Fe] abundance ratio. 

Considering first the primordial stars in the middle and upper panels of Fig. [9l there does 
not seem to be any obvious trend of increasing ^^Mg and ^^Mg fraction relative to total Mg with 
increasing [Fe/H], although it may be that any (small) trend is disguised by the errors in the 
determinations. Indeed the location of the primordial-star points in the lower panel of figure. 



®We assume that the [Fe/H] values are also free of systemat i c effect s. This is supported to some extent by the fact 
the mean iron abundance found for NGC 6752 in lYong et al.l (|2003al ). [Fe/H] = -1.62, is in reasonable accord with 



the mean iron abundance, [Fe/H] = -1.52, of the six NGC 6752 red giants analysed in Norris fc Da Costal (|l995bl ) 
Similarly, the [Fe/H] determinations for M4 in llvans et al.l ( 1999 ) and I Yong et al.l (2008a|) are in excellent agreement 
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[Fe/H] 

Fig. 9. — Relative abundance fractions for ^^Mg (lower panel), ^^Mg (middle panel) and ^^Mg 
(upper panel) compared to the total Mg abundance plotted against iron abundance [Fe/H]. Symbols 
are as for Fig. HI 
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Fig. 10. — Relative abundance fractions for ^^Mg (lower panel), ^^Mg (middle panel) and ^^Mg 
(upper panel) compared to the total Mg abundance plotted against [Al/Fe], the aluminium to iron 
abundance ratio. Symbols are as for Fig. [H 
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where effectively the ^^Mg and ^^Mg determinations are averaged, does support the existence 
of such a trend: making no attempt to distinguish ^^Mg from ^^Mg, the relative isotopic ratios 
appear to change smoothly from ~84:7:7 at [Fe/H] « -1.8 to -^ 70:15:15 at [Fe/H] ss -0.7 dex. Such 
a trend is qualitatively consistent with the nucleosynth etic predictions of inc reasing ^^Mg/^^Mg 
and 26Mg/24Mg with [Fe/H] fe.g. jKobavashi et ahlboilh . although as noted in lYong et all (|2003al . 



20061 ). for example, the isotopic ratios for these primordial globular cluster stars are significantly 



higher t han for field stars of comparable [Fe/H] abundance, despite being otherwise chemically 
similar. lYong et al.l (|2003al . l2006l ) postulate the existence of a generation of relatively massive 
AGB stars in the proto-cluster environment. These stars uniformly enrich the gas, from which 
the current primordial population stars form, in the heavy isotopes of Mg to a level above that 
generated by metal-poor supernovae, which is represented by the Mg isotopic ratios for the general 
halo field population at tha t metallicity . Given the s imilarity of our results for the cj Cen and M4 
primordial stars to those of lYong et al.l (J2003al . l2006l ). it is apparent that this postulated scenario 
could also apply to the proto-oj Cen and the proto-M4. 

Turning now to the full sample of stars, it is noteworthy from the middle panels of Figs. [9] and 
[To] that there does not seem to be any obvious difference between the primordial and the extreme 
stars as regards the ^^Mg fractional abundance, regardless of [Fe/H] and regardless of the [Al/Fe] 
ratio. T he lack of co r relation betw een ^^Mg fractional abundance and [Al/Fe] has already been 
noted by lYong et al.l (J2003al . l2006l ) and it is apparent that the result applies to the to Cen and 
M4 stars observed here as well. Further, the intermediate population stars (predominantly from 
NGC 6752) also have similar values, and a single ^^Mg to total Mg fraction of order 10% appears 
consistent with all the determinations. Indeed for the 34 stars shown in the figures, the mean 
^^Mg to total Mg fraction is 9.8% with a standard deviation of 4.6%, which is comparable to the 
uncertainties in the individual determinations. In other words, in these globular cluster stars, which 
cover approximately 1 dex in [Fe/H] and 1.3 dex in [Al/Fe], and with the current uncertainties, 
the ^^Mg relative abundance appears not to be infiuenced by whether the star has been formed 
from the original cluster material, i.e. the primordial population, or from significantly "polluted" 
material as in the extreme population stars, which are presumed to incorporate material that has 
been processed through the Mg-Al cycle. This lack of any obvious difference between the different 
populations is surprising as, a priori, one would not expect the relative ^^Mg abundance in Mg-Al 
cycle processed material to be es sentially identical to that f or unprocessed material. In this context 
however, it is worth noting that iMelendez &: CohenI (120091 ) do find a small difference in the ^^Mg 
relative abundance between the CN-weak (primordial) and the CN-strong (intermediate) stars in 
the metal-rich globular cluster M71. Specifically, for 5 CN-weak stars, the mean value of the ^^Mg 
to total Mg relative abundance is 5.2 zb 0.4 percent, while for 4 CN-strong stars the mean ^^Mg 
relative abundance is slightly higher, 7.3 it 0.4 percent. In both cases the error given is the standard 
error of the mean. 

The situation for the ^^Mg fractional abundance, however, is markedly different. The upper 
panels (and, given the lack of variation in ^^Mg, the effectively mirror images in the lower pan- 
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els) of Figs. [9] and [TOJ reveal clear dependencies on both [Fe/H] and on [Al/Fe]. In particular, 
Fig. |9] shows that the production of ^^Mg is notably less at higher overall abundance. This is 
consistent with the fact that, in general, the extreme and intermediate population stars in glob- 



ular clusters with larger [Fe/H] values show lower enhancements in 



metal-poor systems (e.g.. iNorris &: Da Costal Il995bl: lYong et al. 



2008a 



Al/Fe1 compared to more 



Melendez fc CohenI l2009l : 



Johnson fc Pilachowski 



201C 



Da Costa &: Marino 



2011 



Carretta et al. 



20131 ). This is usually in- 



terpreted (e.g.. I Johnson fc Pilachowskill2010l ) as a consequence of lower temperatures at the bottom 
of the convective envelope in the more metal-rich relatively massive AGB stars t hat are assumed 



to be the source of the [Al/Fe] enhancements. Indeed, at least as regards u Cen. iD'Antona et al 



(J201ll ) claim that the smaller amou nt of [Al/Fe] enhancernent s een at higher [Fe/H] values is con- 
sistent with the predictions of the I Ventura fc D'Antonal (|2009l ) AGB star models. The reduced 
production of ^^Mg at larger [Fe/H] seen in Fig. [9] is then presumed to be consistent with the same 
interpretation, given the strong connection between the ^^Mg fractional abundance and [Al/Fe] 
seen for the extreme stars in Fig. [TOl 

As regards the data presented in the upper panel of Fig. [TOl and the implicit mirror image in 
the lower panel given the lack of variation in ^^Mg, it appears that there are effectively two regions 
in the relation between ^^Mg and [Al/Fe]. The first is occupied by the primordial and intermediate 
population stars. Here the [Al/Fe] ratio increases by approximately 0.6 dex from the low value in 
the primordial stars but the ^^Mg fraction remains essentially constant, as does the ^^Mg fraction. 
The second region is occupied primarily by the extreme stars though it includes some intermediate 
populations stars. Here we see that the ^^Mg abundance fraction is tightly coupled with that of ^^Al 
as the [Al/Fe] ratio rises from ~0.6 dex to the value of '--'1.2-1.3 in the most extreme population 
stars. At the same time the ^^Mg fraction rises from approximately 10% to almost 40% while the 
^^Mg fraction drops to '^50%. 

The first region of the relation can be understood if it is assumed that the increased Al 
abundance is generated from ^^Mg via the Mg-Al cycle. For example, if we assume the following 
parameters: [Fe/H] = -1.62, [Mg/Fe] = +0.5, [Al/Fe] = 0.0 and Mg isotopic ratios of 80:10:10 (i.e., 
those of the NGC 6752 primordial stars), then, ignoring any dilution, converting just 10% of the 
Mg atoms present to ^ Al results in an increase in [Al/Fe] of 0.7 dex while the (total) [Mg/Fe] 
is reduced by only 0.05 dex. Such a small change in [Mg/Fe] is within the current li mits on the 



variat ion in [Mg/Fe] for the NGC 6752 primordial and intermediate population stars (JYong et al 



2003al ) . Further, assuming that the abundance of the more massive Mg isotopes remains relatively 
unaltered, then the isotopic ratios change only to 78:11:11. Hence in this region of the diagram, 
the increase in [Al/Fe] and the lack of any significant change in the isotope ratios is consistent with 
the conversion of a small amount of ^^Mg to Al in the Mg-Al cycle. We note that the data shown 
in Fig. [To] appear to require the consumption to be primarily of ^^Mg; if the synthesised Al atoms 
in our example come entirely from ^^Mg or ^^Mg then that would produce significant changes in 
the ^^Mg or ^^Mg fractional abundances inconsistent with the observations shown in Fig. [101 



In the second region of Fig. [10] the production of Al and ^^Mg are apparently closely related. 



27 



This is somewhat surprising, as in the conventional treatment of the Mg-Al cycle, ^^Mg is destroyed 
by proton captures (to form ^^Al) and is not expected to be made in any significant amount. 
Specifically, the proton-capture reaction on ^^Mg produces ^^Al, which, in the ground-state, has a 
long /3"''-decay time of ~0.72 Myr. Consequently, the cycle progresses from ^^Al primarily through 
the pathway ^^Al(p,7)^'''Si(/3"'")^^Al and no ^^Mg is synthesised. ^^Al does, however, possess an 
excited state isomer, ^^Al™, for which the /3~''-decay time is short (6.3 sec, cor nparable to the 0~^ - 
decay time of ^^Si) thus providing a channel to produce ^^Mg. Unfortunately, llliadis et al.l ( 201ll ) 
indicate that the rate for the reaction ^^Mg(p,7)^^Al'" is approximately a factor of five lower than 
the rate for the equivalent reaction involving the ground-state of ^^Al at the temperatures likely 
to apply in the hot-bottom burning regions in relatively massive AGB stars. Consequently, the 
generation of ^^Al via the Mg-Al cycle is not expected to show the strong correlation between ^^Al 
and ^^Mg seen in Fig. [TOl 

In the context of the heuristic example outlined above, if we now suppose 50% (rather than 
10%) of the available Mg atoms are processed, including half of the original ^^Mg (but not ^^Mg), 
and that the output of the processing is 80% ^'^Al and 20% ^^Mg (or equivalently 20% ^^Al that is 
assumed to ultimately decay to ^^Mg) then the resulting product has [Mg/Fe] = 0.28, [Al/Fe] = 1.24 
and ^^Mg, ^^Mg, and ^^Mg fractions of the total Mg abundance of 58%, 8% and 33%, respectively. 
These values are consistent with the data of Fig. [10] for the extreme stars. Moreover, the change 
in [Mg/Fe] from +0.5 to +0.28 is co nsistent with the l ower [Mg/Fe] abundance ratios seen in the 



NGC 6752 extreme population stars (JYong et al. 



2003a 



i, see also Table [T]) . While these calculations 
are admittedly ad hoc, they do serve to illustrate that in order to understand the (^^Mg, [Al/Fe]) 
relation revealed by the extreme population stars in Fig. [TOl synthesis of both ^^Al and ^^Mg from 
^"^Mg would seem to be required. 

The discrepancy between the observational results presented in Fig. [10] and cu rrent theory is 



highli ghted by a comparison with the AGB-star evolutionary models presented in IVentura et al. 



(|201ll ). Their "reference model" is for a 6 Mq star with a composition typical of intermediate 
metalhcity globular cluste rs ([a/Fe] = +0.4 a. nd Z = 0.001 or [M/H] « -1.3) evolved through the 



complete AGB phase (see lVentura et al.ll201ll . for details). The temperature at the bottom of the 
convective envelope reaches sufficiently high values that the Mg-Al cycle is initiated and over the 
course of the AGB evolution, i.e., prior to the loss of most of the envelope via mass loss, in the 
model the surface Mg abundance declines by about 0.35 dex while the surface Al abundance rises 
by ~1 dex, which is in approximate agreement with the observations. However, in this model 
26]\/[g ig destroyed throughout the AGB evolution by the p-capture reaction and not synthesised, 
so that there is little real change, from its initial low value, in the ^^Mg fraction of the total 
surface Mg abundance (which decreases during the evolution). In contrast, with the reaction rates 
and the treatment o f thermal pulses, 3rd dredge-up and convective mixing used in the model (see 



Ventura et al. 



2011 



, for details), the total surface Mg abundance becomes dominated by that of 
^^Mg, so that by the end of the evolution the ^^Mg fraction of the total Mg abundance exceeds 

It is very difficult to see how these surface fractional abundances can 



(IVentura et al. 



2011 
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be reconciled with the observations of the extreme population stars shown in Fig. [TUl even allowing 
for dilution of the AGB star surface abundances with pristine material, which is not depleted in 
overall Mg abundance and which has fractional isotopic abundances like those of the primordial 
stars. 

We speculate here on two possibilities that might give rise to a higher ^^Mg fractional abun- 
dance in the most Al-rich stars assuming, given the strong correlation seen in Fig. \lU[ that the 
general contex t is th at of the Mg-Al cycle and AGB star evolutionary calculations such as those of 
Ventura et al.l (J201lh . In the first scenario we postulate that the ^^Mg is sy nthesised as ^^Mg is de- 
pleted. To achieve this we suggest, despite the results of llliadis et al.l (J201ll ). an increase in the rate 
for the reaction ^^Mg(p,7)^^Al™ that produces the excited state isomer ^^Al™. As noted above this 
isomer rapidly /3~''-decays to produce ^^Mg. An increase in this rate and in that for ^^Mg(p,7)^^AF 
wo uld also decrease the ^^Mg fractional abundance reducing the "over-production" of ^^Mg seen in 
the lVentura et al.l ( 201ll ) models. In this scenario the rate for the reaction ^^Mg(p,7)^^Al may also 
have to be reduced to ensure the synthesised ^^Mg is not entirely processed further to ^^Al. 

In the second scenario we postulate that ^^Al (as distinct from ^^Mg) is synthesised as ^^Mg 
is depleted. Here we note that once Mg-Al cycle processing has ceased, material rich in ^^Al will 
eventually show an increase in the ^^Mg relative abundance as the ^^Al radioactively decays to 
^^Mg. To achieve this we suggest a significant increase in the ^^Mg(p,7)^®AF reaction rate (and 
probably also a reduction in the rate of the ^^Al(p,7)^ Si reaction, a lthough ^^Al still nee ds to be 
produced from the /3"''-decay of ^^Si). Here, again in the context of the lVentura et al.l (J201ll ) models. 



the destruction of Mg does not primarily produce ^^Mg but rather instead results in significant 
amounts of ^^Al (as well as ^^Al). This "over-produced" ^^Al is then subsequently observed as 
^^Mg, the radioactive decay product of ^^Al. The suggested rate increase would also serve to 
alleviate the discrepancy between the theoretical predictions and the observations as regards the 
25 Mg relative abundance. 



It is worth noting that the exploratory H-burning nucleosynthesis calculations of lPrantzos et al. 



([20071), which were carried out at a variety of fixed temperatures rather than in the context of a full 
evolutionary stellar model, and which employed the best available reaction rate data, were able to 
reproduce the Mg isotope results for the globular cluster NGC 6752 given in lYong et al.l (J2003al ) as 
well as the distributions of the O, Na and Al abundances. The agreement required quite specific 
conditions in their relatively simple analysis: a temperature of approximately 75MK (below 70MK 
^^Mg is hardly affected while above 80MK it is rapidly destroyed) and for the extreme stars, a 
dilution of the processed material via the incorporation of approximately 30% material having the 
primordial composition. They also point out that the "exces s" of ^^Mg seen in the NGC 6752 stars 
can only be understood if it is originally produced as ^^Al ( Prantzos et al.l 120071 ). Whether these 
fixed temperature results apply to the more complex situation in real stars is a question that is 
beyond the scope of this paper. 



We note only that our results for the Mg isotopic ratios in a; Cen extreme population stars 
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reinforce the discrepancy betw een o bservations a nd ste ll ar evo lutionary theory already apparent in 



the results of IShetrond (jl996l ) and lYong et al.l (|2003al . l2006l ) . A resolution of the discrepancy is 



likely to lead to further insight into the complicated nucleosynthetic processes that occurred early 
in the lifetime of most globular clusters and in that respect further determinations of magnesium 
isotopic ratios in additional globular cluster stars would be very worthwhile. The data presented 
in Figs. [9] and [lO] suggests strongly that the stars to target in such a program, those most likely 
to have significantly non-solar Mg isotopic ratios, are the stars with the highest enhancements of 
[Al/Fe] in globular clusters with intermediate or lower metallicities, i.e. [Fe/H] < -1.4 dex. 
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Table 2. Globular Cluster Red Giant Mg isotope ratios 



Star 


Region 


25Mg/24Mg 


(T 


26Mg/24Mg 


a 


24Mg:25Mg:26Mg 


uo Gen ROA 94 


5134. 6A 


0.000 


0.031 


0.190 


0.023 


84:00:16 




5138. 7A 


0.006 


0.054 


0.029 


0.038 


97:01:03 




5140.2A 


0.139 


0.091 


0.000 


0.066 


88:12:00 




Mean 


0.027 


0.016 


0.106 


0.012 


88:02:09 


uj Gen ROA 43 


5134.6A 


0.161 


0.015 


0.225 


0.010 


72:12:16 




5138.7A 


0.167 


0.021 


0.148 


0.011 


76:13:11 




5140.2A 


0.066 


0.023 


0.194 


0.016 


79:05:15 




Mean 


0.137 


0.006 


0.190 


0.004 


75:10:14 


u Gen ROA 100 


5134.6A 


0.368 


0.068 


1.005 


0.065 


42:16:42 




5138.7A 


0.011 


0.086 


1.087 


0.097 


48:01:52 




5140.2A 


0.497 


0.115 


0.461 


0.072 


51:25:24 




Mean 


0.282 


0.029 


0.835 


0.025 


47:13:39 


UJ Gen ROA 132 


5134.6A 


0.148 


0.012 


0.279 


0.010 


70:10:20 




5138.7A 


0.064 


0.017 


0.204 


0.011 


79:05:16 




5140.2A 


0.176 


0.020 


0.061 


0.012 


81:14:05 




Mean 


0.130 


0.005 


0.188 


0.004 


76:10:14 


UJ Gen ROA 150 


5134.6A 


0.908 


0.091 


2.189 


0.106 


24:22 


53 




5138.7A 


0.124 


0.042 


0.408 


0.031 


65:08 


27 




5140.2A 


0.203 


0.053 


0.234 


0.029 


70:14 


16 




Mean 


0.312 


0.019 


0.550 


0.013 


54:17 


30 


w Gen ROA 201 


5134.6A 


0.163 


0.013 


1.095 


0.036 


44:07:48 




5138.7A 


0.278 


0.025 


0.036 


0.004 


76:21:03 




5140.2A 


0.497 


0.036 


0.000 


0.008 


67:33:00 




Mean 


0.259 


0.007 


0.098 


0.002 


74:19 


07 
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Table 2 — Continued 



Star 


Region 


25Mg/24Mg 


a 


26Mg/24Mg 


a 


24Mg:25Mg:26Mg 


u Cen ROA 248 


5134.6A 


0.165 


0.014 


0.481 


0.017 


61 


10:29 




5138.7A 


0.078 


0.009 


0.082 


0.004 


86 


07:07 




5140.2A 


0.241 


0.027 


0.230 


0.012 


68 


16:16 




Mean 


0.134 


0.005 


0.173 


0.003 


77 


10:13 


M4 L1514 


5134.6A 


0.058 


0.007 


0.226 


0.008 


78:05:18 




5138.7A 


0.111 


0.014 


0.062 


0.010 


85:09:05 




5140.2A 


0.085 


0.013 


0.109 


0.014 


84:07:09 




Mean 


0.078 


0.003 


0.142 


0.003 


82:06:12 


M4 L2406 


5134.6A 


0.054 


0.022 


0.226 


0.015 


78:04:18 




5138.7A 


0.053 


0.038 


0.243 


0.025 


77:04:19 




5140.2A 


0.129 


0.040 


0.073 


0.022 


83:11:06 




Mean 


0.073 


0.010 


0.185 


0.007 


80:06:15 


NGC 6752 mgO 


5134.6A 


0.163 


0.019 


0.495 


0.014 


60:10:30 




5138.7A 


0.029 


0.027 


0.535 


0.019 


64:02:34 




5140.2A 


0.285 


0.038 


0.339 


0.023 


62:18:21 




Mean 


0.148 


0.009 


0.467 


0.006 


62:09:29 



